tridiagonal solver matlab